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We present the calculations of the supersymmetric QCD corrections to the total 
cross sections for single top production at the Fermilab Tevatron and the CERN 
Large Hadron Collider in the minimal supersymmetric standard model. Our results 
show that for the s-channel and t-channel, the supersymmetric QCD corrections 
are at most about 1%, but for the associated production process pp — ► tW, the 
supersymmetric QCD corrections increase the total cross sections significantly, which 
can reach about 6% for most values of the parameters, and the supersymmetric 
QCD corrections should be taken into consideration in the future high precision 
experimental analysis for top physics. 
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I. INTRODUCTION 



The search for single top quark production is one of the major aims of both the Fermilab 
Tevatron and the CERN Large Hadron Collider (LHC) [l], E|, because it can probe the 
electroweak sector of the Standard Model(SM), in contrast with the dominant QCD pair 
production mechanism, and provide a consistency check on the measured parameters of the 



top quark in the QCD pair production process [3|, . Further more, the mass of top 

quark is comparable to the electroweak(EW) symmetry breaking scale, and it can play a role 
of wonderful probe for the EW symmetry breaking mechanism and new physics, especially 
the minimal supersymmetry standard model(MSSM) as a very popular model bey ond the 
SM, via non-standard couplings , [10, [ill, Il2|, Il3j], loop effects 

etc., in the processes for single top quark production. 

At the LHC single top quarks are produced primarily via the t-channel 
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18|, 



Q 



q' + t, 



the quark annihilation process(s-channel) 20l. |21|. 

q + q' -»• t + b, 
and the associated production process 0, Q , 

g + b->t + W~, 

which can reliably be predicted in the SM, and their leading order(LO) results are summa- 
rized in Table [I] [2| . As a high-luminosity LHC would run soon and allow accurate measure- 

TABLE I: The LO results for single top production at the LHC 



process: 


t-channel 


s-channel 


Wt 


aijpb): 


156 ±8 


6.6 ±0.6 


14.0l|* 
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ments of these cross sections with a statistical uncertainty of less than 2% HI, [2 
at this level of experimental accuracy, calculations of the radiative corrections are necessary 
to test the predictions of the SM and to investigate loop effects arising from new physics. 
The QCD corrections to the total cross sections of the three channels are about —10%, 



+50% and +50%, respectively 
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29| . For t- channel, the SM EW corrections 



are about 10% and the supersymmetric(SUSY) EW corrections are a few percent [30J at 
the LHC, for s-channel, the combined effects of SUSY QCD, SUSY EW and the Yukawa 
couplings can exceed 10% for small tan/3(< 2) but are only a few percent for tan /3 > 2 at the 



Tevatron 
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171 ]. and for the associated production at the LHC, the SM EW and the SUSY 



3l|. But at the LHC, the SUSY QCD 



corrections have been calculated and analyzed in Ref. 
corrections to the three channels have not been calculated yet. As the associated production 
process involves the QCD coupling in initial state, obviously, the SUSY QCD corrections are 
significant for this channel. In order to fill in the blanks in the relevant radiative corrections 
to the channels considered here, in this paper we present the calculations of the SUSY QCD 
corrections to the three channels at the CERN LHC, mainly concerning the associated pro- 
duction channel, and also give the updated numerical calculations to the s-channel process 
at the Tevatron 32 ] . 



Our paper is organized as follows. In Section II we will give the analytic results in terms of 
the well-known standard notation of one-loop Feynman integrals for the associated channel, 
s-channel and t-channel. In Section III we will present our numerical results with discussions 
of their implications. 

II. ANALYTIC RESULTS 
A. Associated Production 



In our calculations we use dimensional reduction to control all the ultraviolet divergences 
in the virtual loop corrections. We adopt the on-mass-shell renormalization scheme 33J for 
the top quark mass and the wave functions renormalization, and set all the other quark 
masses as zero. The QCD coupling constant g s is renormalized in the modified minimal 
subtraction scheme (MS) except that the divergences associated with the top quark and 



colored SUSY particle loops are subtracted at zero momentum |34l] . Denoting ip q0 , m qQ , g sQ 
and A^q as bare quark wave functions, quark masses, strong coupling constant and gluon 
wave function respectively, the relevant renormalization constants are then defined as 

m q0 = m q + 5m q , 



i 



A, = {l + 5Z AA fl 2 A^ 

g s0 = (1 + 5Z g )g a , (1) 

After calculating the self-energy and vertex diagrams in FigfTJ we obtain the explicit 
expressions of all necessary renormalization constants as follows: 

a,(M 2 ),A,l , „ JV, , M i, 1 , Ml, 1 , , M L, 

SZ ° = -47 i (f'7-^ + ln < 4 ^ + 3 £ I2 hl <7? )+ E 12 ln ^ )} ' 

u=u,c,t d=d,s,b 

SZ AA = -^S,(0, M|, Mi) + 3 ^B m (0, Ml M?) 

E E *£. M D-^tt «|. 

Mi d i 

2 2 

i i 

5Z U = - 2 ^t B ^ M lK)uk^ 



3?r Z r ,-^~t,'-g>'-t i >''WW) 3yr Z_^ p 2 



2a s m 2 v 9 



3/T ^— ' dp 2 

M q { l R) = -^J2 b ^ M l M l) u h» u lm 9 = M, «, c, 6, (2) 

i 

where /?o = [— \{N + 1) — |(n/ + 1)], -B^... are the two-point functions, are the mixing 
matrices of the squarks, 7^ is the Euler constant, and /1 is the renormalization scale. 
Including the SUSY QCD corrections, the renormalized amplitudes can be written as 

M r l = Mq + Mt + Mi unt , (3) 

where Mq is the LO amplitude, M^ ir represents the SUSY QCD corrected amplitude from 
the one-loop self-energy, vertex, and box diagrams, and M^ ount is the corresponding coun- 
terterm for the self-energy corrections and vertex corrections, respectively. 

The related Feynman diagrams which contribute to the LO amplitude are shown in 
FigfS] (a) and (b). The LO amplitude M$ is given by 



Mq = Mf + M i 

= -^-[ 1 -(2A 10 -A 7 -2A 8 ) + -^ 
V 2 sin Q w s t-m t 



[-(2A W -A 7 - 2A 8 ) + ——^(2A 15 - 2A 9 - A 7 )], (4) 
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where are the elements of CKM matrix, s, t, u are the Mandelstam invariants, which are 
defined as, 

s = (ki + k 2 ) 2 = (k 3 + h) 2 , 
t = (h- hf = (k 2 - k A )\ 

u = (h - h) 2 = (k 2 - k 3 ) 2 , (5) 

where k\ and k 2 denote the momentum of the incoming particles, k 3 and k± the outgoing 
particles, while A m are the reduced standard matrix elements given by 



Al,22 — 


(4 


■el)u h (k 3 )P R , L u c (k 2 )(T a ) bc , 


A 2 , 3 = 


(3 


■ k 2 )(el k h2 )u b (k 3 )P R u c (k 2 )(T a ) bc , 


^4,5 = 




■ h)(el k 2:1 )u b (k 3 )P R u c (k 2 )(T a ) bc , 


^6,17 = 


u b (k 3 )P R , L ft ft 4 u c (k 2 )(T a ) bc , 


Af,27 = 


u\h)P R)L ft ft An c (A; 2 )(T a ) 6c , 


^8,9 = 


{el 


■k 2 , 3 )u b (k 3 )P R ft A u c {k 2 )(T a ) bc , 


^10,32 = 


(c? 


■et)u b (k 3 )P R , L Mu c (k 2 )(T a ) bc , 


4u,12 = 


{<• 


' k 2 )(el k 2 , 3 )u b (h)P R fau c (k 2 )(T a 


^-13,14 = 


«■ 


' fci)(e?- k 2>3 )u b (k 3 )P R fau c (k 2 )(T a 


^4-15,16 = 


(el 


K 2 )u b (k 3 )P R . ftu c (k 2 )(T a ) bc , 


^18,19 = 


(el 


h, 2 )u b (k 3 )P R ft fau c (k 2 )(T a ) bC) 


^20,21 = 


(el 


■k 2 )(el-k h2 )u b {k 3 )P L u%k 2 )(T a ) bci 


^23,24 = 


(el 


■k 3 ){el-k 2 ^u b {k 3 )P L u%k 2 ){T a ) bc , 


^25,26 = 


(el 


K 2 )u b {k 3 )P L ftu c (k 2 )(T a ) bc , 


^28,29 = 


(el 


K 2 )u b {k 3 )P L ft Mu c (k 2 )(T a ) bc , 


^30,31 = 


(el 


■k 2 , 3 )u b (h)P L ftu c (k 2 )(T a ) bc , 


^33,35 = 


(el 


k 2 ,i)(elk 2 )u b (k 3 )P L / k 1 u c (k 2 )(T a ] 


^34,36 = 


(el 


k 3 )(el-k 2tl )u b {k 3 )P L Mu c (k 2 )(T a [ 


^37,38 = 


(el 


■ k 3 , 2 )u b (k 3 )P L ft Mu c (k 2 )(T a ) bc . 



The relevant Feynman diagrams of the SUSY QCD corrected amplitude M^ ir are shown 
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in Figj3j and M^ ir can be written as 

Mt = Mi lf + Mttex + M L, (7) 

where M^j, M^ ertex and M^ ox come from self-energy diagrams, vertex diagrams and box- 
diagrams as shown in FigJHl respectively. Their explicit expressions are given by 

38 38 38 

M^elf = ^2 fm*A m , M^ ertex = fm^m, ^box = ^ fm^m, (8) 

m=l m=l m=l 

where f^ 1 -^, f^, are the form factors, which are given explicitly in Appendix. 

The counterterms M^ unt , the corresponding diagrams of which are shown in Figj4j can 
be written as follows 

MLnt = 5Mf elf + 5M* rtex , (9) 

with 

5Mf elf = 5M] elf + 5M 2 self 

= 6 " sVtb a { -^[m t (SZ* - 5Zl) + 2Sm t ] - ^l(A 7 - 2A S + 2A W ) 
2\j2sm6 w t - mf s 

+ 7T 2 2 , 2 (A 7 + 2A 9 - 2A 15 )[(m 2 t - t)SZf + 2m t 5m t ] }, (10) 
[t m t ) 

4 

$ M vertex = ^2 5 M vertex, ( U ) 
n=l 

SMLtex = ~ T V !\ (Ar + 2 A 8 - 2A 1Q ) (6Z» + 6Z») , 
2\J 2s sin v w 

SMLtex = ~ 7E . e9 : V u ^(A 7 + 2A 9 -2A 15 )(5Z b * + 5Z t t), 

2a/2 sm 9 w [t - mf) 

SMLtex = ~ 7= . ^fl" V A2m t A 17 {5Zt-5Zt) 



2\/2sin 6 w {t - 



-(A 7 + 2A 9 - 2A 15 )(25Z g + 25 Z% + bZ AA )\ 



SMLtex = ~ e J sV !\ (2SZ g + 5Z AA + 25Z b L b )(A 7 + 2A 8 - 2A W ). (12) 
2 V 2s sm Uw 



The partonic cross section can be written as 

•i 1 

a 



l d ^ AV2 i^-i 2 = r d w |M " j2 ' (i3 » 

where A = (s - m\ + m 2 ^) 2 - 4sm 2 y , £± = |[m 2 + m 2 ^ — s ± A 1/2 ], and \M^ en \ 2 is the 
renormalized amplitude squared given by 

JmJJ = 1 Mi 1 2 + 2Re^M A [M*. + M? ount ] \ (14) 
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where the colors and spins of the outgoing particles have been summed over, and the colors 
and spins of incoming ones have been averaged over. 

The total cross section at the LHC is obtained by convoluting the partonic cross section 
with the parton distribution functions (PDFs) G g ^/ p in the proton: 

a = / dxi / dx2[G g / p (xi,/if)G b / p (x 2 ,iJif) + (x x x 2 )]a(rS), (15) 

Jto Jto/x\ 

where /!/ is the factorization scale and S = (P\ + P2) 2 , Pi and P 2 are the four-momentum of 
the incident hadrons, r = ( m ^+ mt ) ; T = xix 2 , and x±, x 2 are the longitudinal momentum 
fractions of initial partons in the hadrons. 



B. s- channel and t- channel 



For convenience, we first define the reduced standard matrix elements Fj as follows: 

F 1>2 = v(k 2 )P RjL ^u(k 1 )u(k 3 )P R -f IJ ,v(k4), 

F 3 ,4 = v(k 2 )P R , L u(ki)u(kz)P R fciv(k 4 ), 

^5,6 = u(k 3 )P RtL v(k 4 )v(k 2 )P R faufti), 

F 7 , 8 = v{k 2 )P R>L u{ki)u(k 3 )P L v{k4), 

F 9 = vik^PR^uikjuik^PLj^ih), (16) 

which will appear in the amplitudes of s-channel and t-channel below. 

For s-channel, the diagrams which contribute to the LO amplitude Mq are shown in FigfS] 
(c). The LO amplitude Mq is 

M o = - £ • 2n f ^ 2 A (17) 
sm 9 w (s - m z w ) 



q—u,c 
q 1 — d,s,b 



The virtual corrections M^ ir contains the radiative corrections from the one-loop vertex 

diagrams, which are shown in Figj5] (a) and (b), and we can write M% ir as: 

9 

M s mr = £ &F m , (18) 

m=l 

where /4( m = 1; 2, ■ • • ,9) are form factors, which are given explicitly in Appendix. 
The corresponding counterterm can be written as 

M'eaunt = ~ £ . ™^ 2 A 52!» + 52% + 52% + 52$')^, (19) 



q — u,c 
q 1 — d,s,b 



s 



where the expressions of 5Z™(i = u, d, s, c, b, t) are shown in Eq.([2]). 

According to the crossing symmetry, we have the similar expressions in t-channel as in 
s-channel. We replace the variable s by t in Eq. ffTTl) . (fT8l) and (fT9l) . then use the different 
summation over quark flavors, and change the indices of the quarks in the initial and final 
states. For example, 

9 

q,q' {qq>} m=l q,q> { qq >} 

M 'count = M count (s — > t, 2^ ~~ > 

q,q' {qq'} 

= - E • Z a Y qq \^ z i q + sz i q ' + 5Z ? + 6Z *w (2°) 

*—f sin &w{t — m - 



{qq'} 



Wl 



with 



F' m = F m {v{k 4 ) -> u(k 2 ),v(k 2 ) -> u(k 3 )). 

Here the index pair {qq'} takes on the flavors {ud, us, ub, cd, cs, cb}. The other formulas for 
cross sections are the same as in associated production. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we present the numerical results for the SUSY QCD corrections to the 



three channels of single top production at the 
numerical results for s-channel at the Tevatron 



HC. For comparison, we also present the 
321. In our numerical calculations, we use 



the following set of the SM parameters 35J: 



m t = 175GeV, a ew (M w ) = 1/128, a s {M z ) = 0.118, 

and all light quark masses are set to be zero, and the CKM matrix elements are taken to be 
the values shown in Ref. 35]. 



The running QCD coupling a s (Q) is evaluated at two-loop order 36J, and the CTEQ6M 



PDFs 



371 | are used throughout this paper to calculate cross sections. For simplicity, we 
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neglect the b-quark mass. We choose /x r = fi f = m t + m\y for the renormalization and fac- 
torization scales in associated production and choose fi r = fif = m t for the renormalization 
and factorization scales in the other two channels. 

Besides, the values of the MSSM parameters taken in our numerical calculations are 
constrained within the minimal supergravity scenario(mSUGRA) 38|, in which there are 
only five free input parameters at the grand unification where M1/2, M , A , tan/3 and the 
sign of /i, where M1/2, M Q , A are, respectively, the universal gaugino mass, scalar mass, and 
the trilinear soft breaking parameter in the superpotential. Given these parameters, all the 
MSSM parameters at the weak scale are determined in the mSUGRA scenario by using the 
program package SUSPECT 2.3 [39J, where we set A = — 200GeV and ji > 0. 



A. Associated Production 



We define the K factor as the ratio of the SUSY QCD corrected cross sections to LO total 
cross sections, calculated using the CTEQ6M PDFs. Figj6] shows the K factors as functions 
of Mg (M1/2) for the associated production process pp — > tW at the LHC for tan/3 = 5, 20 
and 35, respectively. From Fig® we can see that the differences among the results are small 
for different tan/3, and K factors increase with the increasing M g for small Mg{< 160GeV), 
while decrease with the increasing Mg for large M g (> 160GeV), and, in general, the K 
factors are about 1.06. 

In FigJTJ we show the dependence of the K factors on Mi (M ) for different tan f3. FigJTJ 
shows that the K factors have the similar behaviors as those shown in Figj6j and are also 
about 1.06 in general. 

To compare with FigJTJ we present Fig{10] which takes the similar parameters as those 
used in FigJT] but M1/2 = 70GeV, where the gluino mass M g lies in the range 220GeV < 
Mg < 250GeV for all values of Mo and tan /3 we assumed here. FigJTO] shows that SUSY 
QCD corrections are not sensitive to tan/3, which is consistent with FigfTJ 

FigJH] gives the K factors as functions of Mg(Mi/ 2 ) for different M , assuming tan/3 = 
5. In FigJH] we can see that there are large differences between different M when Mg 
becomes small, but these curves approach each other when M g becomes large because of the 
decoupling of heavy gluino(M<j > 450GeV). The K factors are about 1.06 for M g < 500GeV, 
and become small slowly with increasing Mg, but the K factors decrease rapidly when 
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Mg < 150GeV for M = 150GeV. 

In Figj9] we present the K factors as functions of Mf i (M ), assuming tan/5 = 5, and 
M1/2 = 40, 70 and lOOGeV, respectively. From FigJHl we find the similar results as those 
shown in FiglHl i.e. the K factors are about 1.06 for most values of M~ tl considered, and 
become small slowly with increasing M^. 

In FigJTT]we present the LO and the SUSY QCD corrected cross sections as functions of 
renormalization and factorization scales \xj\x^\Xf = \i r = /i, fi = m t + m w ), respectively, 
assuming tan/3 = 5, M = 150GeV, M 1/2 = 70GeV, A = -200GeV and \i > 0. This figure 
shows that the scale dependence of the SUSY QCD corrected total cross section is a little 
larger than that of the LO cross section because of the running effects of the extra a s in 
SUSY QCD corrections. We can recover the LO results of scale dependence by dividing a s 



in the SUSY QCD corrections. After comparison with the NLO QCD corrections 28], we 
can see that if the NLO QCD corrections are also included, 0(a s ) corrections still improve 
the scale dependence. 

B. s-channel and t-channel 

For the s-channel process pp — > tb, in Fig {12] (a)-(b), we display the K factors as func- 
tions of Mg(Mi/ 2 ) and M^(M ), respectively, assuming M = 150,300GeV and Mi/ 2 = 
40, 70, lOOGeV. FigJT2l shows that the K factors are about 1.01 for some favorable parame- 
ters, otherwise, the K factors approach the unit value. 

In Fig{13] (a)-(b), the K factors are plotted as functions of M g (Mi/ 2 ) and M^Mq), re- 
spectively, assuming tan/3 = 5,20,35. FigJTBl shows that the K factors can reach 1.01 for 
light M g and , respectively, and are not sensitive to tan (5. 

In FigHH we display the K factors as functions of M^ for the s-channel process pp — > tb 
at the Tevatron Run II. We find that the K factors can reach 1.01 for small M^, and 
decrease rapidly as Mj increases. These results are consistent with that shown in Fig. 3 of 



Ref . |32| , where the relevant parameters assumed are the same as ones used in our numerical 
calculations. 

For the t-channel process pp — > qt at the LHC, Fig fl~5l shows the K factors as functions of 
M g for tan (3 = 5 and M = 150, 300GeV, respectively. From FigH5l we find that the SUSY 
QCD corrections decrease the total cross sections and the K factors approach the unit value 
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for all parameters assumed here, which means that SUSY QCD corrections are negligible. 

In Fig{16] and Fig(T7| we present the LO and the SUSY QCD corrected cross sections as 
functions of renormalization and factorization scales n/no(fif = \x T = /i, fiQ = m t ) for both 
s-channel and t-channel, respectively, assuming tan/3 = 5, Mq = 150GeV, M1/2 = 70GeV, 
Aq = — 200GeV and fi > 0. Since the SUSY QCD corrections to the two channels are very 
small, they do not affect the scale dependence of the LO results obviously. 



IV. CONCLUSION 



In conclusion, we have calculated the SUSY QCD corrections to the total cross sections 
for single top production at the Tevatron and the LHC in the MSSM. Our results show that 
for the s-channel and t-channel, the SUSY QCD corrections are at most about 1%, but for 
the associated production process pp — > tW, the SUSY QCD corrections increase the total 
cross sections significantly, which can reach about 6% for most values of the parameters, and 
the SUSY QCD corrections should be taken into consideration in the future high precision 
experimental analysis for top physics at the LHC. 
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Appendix 

In this Appendix, we will list the explicit expressions of the non-zero form factors of 
associated production and s-channel. For simplicity, we first define the abbreviations for 



Passarino-Veltman functions 42j below. 

B s = B (s,M 2 g ,Ml), B{ = B±(s, M~, M|), B" = B (s,M 2 g ,M l 

Bl = B (t, Ml Af?), B\ = B^t, Mf , M|), 
Ct 3k ... = C ijk ...(ml s, m 2 w , Ml M? g , M| ), 
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C» jk ... = C ij UOArnlr,Ml,M?,Mi.), 
^. fc ... = ^...K 2 ,t,0,M|,M|,M|), 
C^... = ^ fc ...K,t,0,M|,M|,M|), 
C^... = a ife ...(m t 2 ,t,0,M|,M|,M|), 

(%,.. r, ; ,..(o^.()..\/ ; ;..u 2 ..\/|). 

C* fc ... = C^(0,« l m? > JW? > i^,JW?), 

qJ fc ... = c« fc .>(o,a ) o ) ^ ) ^ ) ^) ) 

C* k ..=C ijk ...(mls^MlMl,Ml), 



D a ijk ... = D ijk ...(0, m 2 w , ml 0, t, s, Mf , Afg , Af| , Aff 



Dj,- fc ... = A ife ...(0, s, m 2 w , u, 0, m 2 , A/g, Aff, Af|, M|), 



For s- channel, we have 



2 8aa s V tb V:„, 



g t^ c r^Tj -3 sm V w {s - m w ) 



q' — d,s ,b 

2 



q' = d,s,b ' J 

fi = - E E ,y ^^(qf + + ctf ), 

q 1 =d,s,b 

ti = - E E o • 2 S M S UIU<U£U£(C* + ^ + Cff), 

g' —a,s,o 

^ 2 8aa a K&K*„/ 



1B1 

q* = d,s ,b 
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/.* = -£ E i ^m^ic? + eg + eg) 

„^r„ ttt, o sin % [s — m w 



q=u,c i =1 
q' — d,s,b 



7 W {a - my,; 

rt TTb TTt*Tjb* ( s~iB i /^y_B , /~<B\ 



+M- g ul 1 u b ]l ul;u^{c» + cf + c*)], 

£rl ^— 1 3 snr 9 w (s- mfy) 



q = u,c 
q f — d,s,b 



»J=1 



^— ' * ciri V(s — m w) 



=1 3 sin 2 ewis-mt 

q f = d,s,b 



8aa s V tb V* q , nBrjt . . 



For associated production, we calculate the individual diagram separately by different 
types, as shown in Figj3]and Eq. (JSJ). The form factors of the self-energy diagrams are 

2 



f seif 1\pla s g s eV th ^ b bit 



2 

3vr singly E^~"~^ (t - m 2 t f 

+(tU^ 1 UH + m 2 t U! 2 U^)Bl}}, 
fS eif _ 2V2a s g s eV tb A b , 



2 

37rs sin 9w 

2 



On Sill t/p^/ [T ''vf ) 



* ' i 

To avoid the very long expressions, we define 

7 

£1) — \ A fUQ 

Jm / J m i 

a=l 

as shown in Figj3j where 

ft = ft = ~ E E 3 L sinV + ^ + ^ " W^M, 

i 3 
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rvl rvl 

J 20 — J 21 



rvl _ 1 rvl _ _ 1 rvl _ _ V" A ^ a ^9sVtbUllU b j l TTh t 

h ~ 2 h ~ 2 /l0 ~ ^2_. 3 S7rS m9 w u H u i^oo, 

* 3 

rvi _ f vi _ _ v- v- 2V2a s eg s V tb Ut 1 U$ TTb 
/i5 - /i6 - 2^2^ 3 7rsin^ ^i^i °i2' 

* i 

fis = fi9 = ~ E E o fl J [M s V$JJ$Ci - m t U 2 U%(C? + - CS, , 

z — ' z — ' ds7rsinc7vv 

» j 

-EE 5 • fl - ["H^i^T(C? + Cn + ^12) - ^^M, 

OSTC Sill t/"^ 

* i 

rvl rvl 

/25 " /26 " 2^ 2^ 37r sin g ^i2 ^12> 

,„i _ I ™i _ _I f ,i _ _ V" V- 2V2a s eg s V tb U t ll U b j l TTh 

/27 - 2 /30 - 2 /32 - 2^2^ 3s7TSm^ ^2^00, 

* J 

/£ = f29 = - E E o fl J [MsUfjJ^Ct - mtU b 3l UH{Ct + C?! + C? 2 )], 

* — ' * — ' OS7T Sill C7W 



2 2 



00; 



/ — 2 2 

rv2 _ 1 rv2 _ _ 1 rv2 _ _ 2\/2a s eg s V tb ^ jjt n b *n b n u r ] 
h ~ 2 h ~ 2 Il5 ~ 37rsin^(t-m 2 )^Y ' 

f v2 _ ^V2a s eg s V tb s—^ y—^ t TT b* TT b TT t* n b 
/l6 " Svrsin^ 2^ 2^ ^i^i ^1^1^12, 

* j 

™2 _ _ 9 ™2 _ _ yu2 _ 4-\/2a? g eg i s Vffc ^-v ^-v . TTb*Tjb rjt*Hb 
/25 " ~ hi ~ 37rsin% ( t _ m 2) U H U H ^2^00, 

2 _ 2\/2a: s eg s Vf6 v-^ rrt rr&*rr& rr**r<& 
i j 

f?+f? = -kf$+f$) = W7^¥t + 5tD( 18C S. +2C S)- 9B S-9^?^)^ 

-m^M^U^Cl + 2M § U t t2 U t l *C c + 9m t U* 2 U**Cl)}, 
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^ 3 + = To m E[( 18m ?^i + 18m * 2c u + 36C7 oo " 18BS - 18tC£ + 4CJ 

12v2vr sin 6^ (t — ) ~^ 

+2m?C* + 2m t 2 Cf 1 - 2iCf 2 + m?C?)E£C/£ + 2m?(9C 1 c 1 + C? + C*)^ t£ 
+2m t M~ g U^U!i(9C e 1 - Cf) + 2M- g m t Uf 1 U!;(9C e 1 - C*)\, 

At + At = ~ ^ iT EW^(9BS - 18C C - 2C d + 9m?CJ - 9tC x e ) 

12v27T sm Vw{t — mf) 

+9M~ g U t ll U t i ;(m* - t)C c + m^U^lSC^ - 95< + 2C d )], 

/ 2 t + fH = -2(f$ + At) = - 7= V * i7 £[M s E&0£(C? - 9C C - 9C0 

3 v 2vr sm 0^ [t — mf) 

+m t Ut 2 Ui*{9C c 12 + Cf 2 ) - mtUlxU^^Ctx + 9C X C 2 + Cf 2 + C? + C d u + 9Cf)], 

zt + a 6 = -2(a 5 + a 6 ) = - *T 9sVtb Q E Mj^ci 

A 5 + A 6 = -^(At + Ao 6 ) = - ^T 9 ^, E " 2C& - 18C&), 
A 5 + A 6 = -TTTp^V XX^Wo - 2C& - 18C&) + rf&l#(9C& + Cf 2 )], 

DV 27TS Sill 

At + At = - - a f 9sVt \ E - 9 ^ f - 9 ^o ), 

DV 27TS sm Vw 

fx 7 = -T^V E E U^UUM^ + m.DSCj), 

At = - - % eg . y " E E + 

6V27rsm% i ~^ 

As before, we define 

3 

f b = 

im / J m i 

0=1 

as shown in FigJ21 where 

2 2 

A 61 = E E U^xUUM.U^ + m^Z% 2 ), 

V27rsm% ij 

A W = ~ E E U^xUUM^Dt, + m^I^) 

V 27r sm W 

IT j j 
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E E UbUSWiVbUTdm + + ^2) + mtU^D^ + D« 2 + L>? 22 )], 
V27T sm 0^ V i 

4?^r E E ^ [M 5 t4^(^ 2 + D« + -D22) + m t Ut 2 U% {D{ 22 + L»« 2 + D a 222 )} , 
V27T sin%Y j 

/f = ~ E E U^UUM.U^ + m^D^), 

V 27r sm 9 W / — 1 *-r 

2v / 27rsm6' H / ^ ^ 

.bi _ 3a s e(7sVf6 TT t T r b* Tr b rr t* n a 

V27T sm 

V27T S1I1 0w ' ' 

/S = _ £ £ f/^M^oV + 

V 27r sm 9 W 

i> 3 

t = - ^ SegsV ' b E E U^U^l(Dl 12 + 0» u + D° 2 + Dl 22 + D° 13 + Z^a), 
V 27r sin^Yy 

" ^T^T" E E ^ ^ ^ + ^122 + ^ + + ^ 22 + + D« 3 + ^3 ) , 

V27T sm 0^ ^ ^ 

m = 3 ^ egsV * E E U^U^l(Dl 33 + D ?23 ), 
V27T sm W 

2 2 

)S = *ro egSV « E E ^K^23 + ^23 + 

/S = - ^T 9 ^ E E ^i^^^(2^oo3 -CI- 21% + m f 2 ^ 3 - ^ 3 ), 
2V27T sm ^77 

~ ^ ^ E E (D« 12 +D« 2 +D° 2 ) -C^D^+D^+D^)}, 

2v2irsm9 w 

13 

2 2 

3 ?^^ E E ^IM^^+D^^ 
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fll = ~ ^ SG9sV n E E U^U^rrnU^ + M~ 9 U? 2 Dt z ), 

V27T SHX&W i j 
2 2 

= 3 * se9sV * b E E C/JiC^r K^?i ^Pn2 + ^12 + ^22) + M s 0* tfg (-Di + D° u + 
V27r sin^Yy 

2 2 

V27T sin^Y^ 

2 2 

= 3 ^ e ^ E E u h u n + ^ + ^ 2 a 22 ) + MfiUft u% (D a 12 + D 2 a + Dy] , 

V 27r sm O w 

fll = - 3 ™ se9sV * E E ^^^(-^1^223 + m^ 3 ), 

V 27r sm 6 W . . 

fll = - J%' e9 : V * E E U^U^mlD^ - W« 3 - A 25 - 2^ + 2%), 
2V27rsm% ij 

_ J^9sV^_ £ ^ u* u*ul U%[{t-ml) (^ 2 +D 2 a +^ 2 )-^-2(^ +D V+^ 02 )], 
2^2^3111^^^ 

/iv 1 = - 9 3 £^ E E ^M^, 



= - 3 ^-^ E E + ^ + Dh) + M~ g UlUf 2 (D« + D« + D})], 



.61 _ 3a s e(7gVffc t TT b* TT b TT t* n a 

V27T sm^ 



/£ = f eg : v * E E ^M^v 

fll = ~ ^o^a E E U^U^D^ + Z% 3 ), 

fll = - ^ seg : v * E E u^u^;(di 12 + 0? u + d? 2 + d? 22 + + ^ 23 ), 

= " £ E E U n U % U i2 U %Wi + D i22 + D a 2 + 2/^ + D a 222 + L>* 23 + D 23 + D<£ 23 ) , 

V 27r sm 0^ 



v / 27rsm6 l H / 

w * 3 

fll = ^ segsVt " E E U^U^(D° 23 + ^ 3 + J}*,,), 
V27T sin # w i ^ 

if = - ^'.^ E E ^ W&?JM, + m^t^oa), 
3v27rsm^ i 

ff = - "f 9 ^* E E ^^[^^(^1 + ^12 + ^23 + ^ + ^3 + D\,) 

3V27T sin^YY 

+m^^Pm + D b 13 + D b 233 + D b 33 + D b 333 + D b 133 )}, 

ff = - E E ^[^^i + -Di 3 ) + m^^(^ 13 + ^33)], 

3V27rsm^w 

/f = E E U^(M^lD\ 3 + mtU*^^), 

3V27rsm9 w ij 

3V27rsm^ ~ j 

/I 2 = E E qjftulv&i&i + ^oos), 

3V2tt sin^Yy 

3v / 27rsm6 l W / ^ ^ 

fll = : E E ^M^So* 

3V27rsm^ 

E E + ^123), 

3V27T sm^Yy 

6 2 = aeg sVtb * * 

3V2irsm6 w Jl tl ll 31 123 

2 2 

/g = - *£ 9aV \ E E ^M^i*PL 3 + D b 23 + D b 33 + 14 + D» 22 + ^ 23 ), 

^f^T- E E ^^^23 + ^3 + ^33), 

3V27rsm% ij 



,62 _ 

ill — 
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3V27T sin^Yy 

/f>2 _ ®s^hVti) ST^ ST^ jjt TTb* TT b TT t* n b 

/ 2 6 2 = - Q ^.^ E E ^*K^(^23 + ^3 + ^33 + + D 

3V27T sm^YY 

+M § t/Af/* 2 *(^ + D\ 2 + D 2 6 3 + ^ + ^ 3 + D b 13 )}, 

2 2 

^ = q T 9 — Q E E UhU^mtU^KDl, + D? 33 ) + M-^U^ + D» 3 )] 
3V27T sin^Yj 



333 + -^133) 



/2I = - a f 9sVtb a E E ^* K^£^o 3 + M~ g UlU%D^ 
3v27rsm^w ■ 

2 2 

ffe2 _ 



f£ = E E Um{m t U^U}lD\ m + M^C/j*^), 

3V2tt sm^YY 

3\/27r sin^ ij 



2 2 

b2 _ a s eg s V tb t b b t b b b 

J25 — „ n: . n 7, / , U jl u il u i2 u j2\ jy 00 + -^002 + -^00: 

2 2 

f b2 _ _ a s^g s vjb ST^ ST^ TTt TTb*Tjb TT t* n 6 
3V27T sm % Y Y 

2 2 

/3 ° " 3V27T sin (V V T 41 

/a? = - a ° egsVt \ E E ^M^03, 

3V2-Ksm9 w 

\r 22 

h2 - 3V27rsin^ 1 1 ' ° 02 ' 

« = o a £ 9aV * E E ^ + 
/ 3 6 I = - t^^ E E ^M^i 6 23 , 

3V27T sm^YY 
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fb2 _ 
J35 ~ 



- *f 9sVtb a E E U^Ut 2 U%{D\ 2 + D\ 22 + ^ 23 
3 V 27r sm W i 



-^223 + -^23 + -^233 + -^123)) 



rb2 _ 
J36 — 



- E E ^M^Pk + d> 3 + ^ 33 ), 

3V2tt sm^Yy 



-fb3 _ 



- E E WKW (^oo + ^001 + ^03) - M § U^*D c J, 

3V2n sm 6 W i j 

-TT^V E E ULU%[M~ g U» 2 U%(D° 23 + D{ 3 + D 3 C + d 3 c 3 : 
3V27T sm 9 W ^ ^ 



m t U b 2 U%(D c n3 + 2D^ 3 + D\ 23 + 21% + D° 3 + D 2 C 33 + D c 3 + 2D 3 C 3 + D 333 )], 
22 



f63 _ 
J3 — 



a 



E E UtiU»l[MsU b j2 UttD{ 3 - m t U b j2 Ut*{D{ 13 + D{ 33 + £>«,)], 

l7r sm 9 W 



3\/2tt sin #n/ 



9 9 

E E ULUfArntU^mi + + 2D\ 13 
3v2tt sm 9w i j 

-M~ g U* 2 UZ(Dl 1 + D c 13 )], 

V 2 2 

- T 9 : - a E E UtxU^mtU^Ul^Dt, + D C 1U + 3D c ll3 + 3D{ 33 + 2D\ 23 + AD 
3v 27r sm 0jy i j 

+D23 + D c 233 + D C 3 + 2D C 33 + D c 333 + D\ + D c 112 + D\ 2 ) 
-M~ g U b 32 Utt{Dl + D\ x + D\ 2 + D- 23 + 21% + D 3 C + D c 33 )], 



rb3 _ 

J A — 



2D{ 13 + D{ 33 + D 



13) 



c 
13 



fb3 _ 
J8 — 



f b3 
J9 



3 v 2ix sm 9 W 

2 2 

- % egsVtb Q E E ^i^r^^(^ c oi - ^oo 3 ), 

3y/2iTsm9 w ^ . 



2 2 

rb3 _ _ a s^g s Vtb rr t Tjb*Tjb jjt* n c 

J10 — Q 7k - q / j / , u il u jl u jl u 11^0021 
3V 2tt sm 9 W ^— ' 



i j 
2 2 



f b3 _ 
ill _ 



f b3 _ 
Jl2 — 



a ' e9 ' v * e E uku b iu b u^Dt 23 , 

3^/2-Ksm9 w L - jL - J 3 J 

13 

3V27T sm Owij 



fS = q a r" 9sVt \ E E v\ x upJ)jJ*{Di 23 

3V27rsm^ i j 



D c 223 + D\ 3 + D 



233/' 
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fb3 _ a * 



til = - a ' e9 ' V * E E UlU^UttiDlu + + ^22 + 2i^ 23 + £ 2 C 23 + ^23 + ^233), 
3V27rsm^ Y j 

2 2 

E E U * U n U k U % ( D oo + 2^o C oi + £002 + ^03), 
2 2 

= - q % egaV * E E ^1 iM^utm* + ^23 + di + ^ 3 ) 

3V27T sm^w Y j 

—mtU^-JJll {D\ X3 + 2DJ 3 + D£a + 2L^ 33 + D c 23 + D 2 C 33 + D c 3 + 2D C 33 + D 3 C 33 )], 

til = ~ - E E ^[Ms^lJ^DU - m t U b n UH{Dl l3 + DJ 3 + ^ 33 )], 

3V27T sin^Yy 

2 2 

/22 s = ~ ^ * E E ^i^rh^x^^So + ^001 + D c m3 ) - M § U^D c 00 ], 

2 2 

= - *£ 9SV \ E E ^^K^i^Pn + + 2D? 13 + ^3 + ^33) 
3v 27T sin 9 W i 

-M- g u b n ut;{Di l + Di 3 )i 

fll = -t^^EE^Xi>*^^ 

3V27rsm^ ^ j 

+ 22^ + 3^33 + + ^33 + D° + 2^3 + ^33) 

-M~ 9 U^(D c n + D{ 2 + D° + 2D\ 3 + L> 2 C 3 + D C 3 + D c 33 )\, 
2 2 

Q a eff a ^& r T t Tjb*Tjb Tjt*( D c , D c , D c _i_ j~) c \ 

J25 — „ /77 . „ / j / j u il u jl u j2 u i2 1-^00 + -^002 T -^003 i~ -^qqiJ, 

3V27rsin^ i j 

2 2 

fb3 _ _ a sZQs Vjb ST^ ST^ jjt Tjb*Tjb TT t* n c 
J2&— „ /k . n ?, 7, u il u jl u j2 u i2 -^001) 

3V27T sm^YT 



t/ 2 2 
f b3 _ a seg s vtb sr^ 77* TT b *TT b tt u n c 

J30 — /7T . ^ u il u jl u j2 u 12^0031 

3v27rsin0w 

/£ = E E ^r^re* + ^ 

m = - ,T sVt l E E umiu%ut*Di Q2 , 

3V27T sin^YY 

/3 = , T sV l E E ^^^^^3- 

3v27rsm^vK 
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2 2 

3V 27r sm 0w i 

2 2 

/sl = - T 9sVtb E E ^^^2^2(^23 + D° n + ^33 + ^ 23 ), 

3V27rsm0 w i y 

/36 = ~ - E E Ut *U%U> 2 U%(D* 223 + D 2 C 3 + -D233 + £>H2 + ^2 + ^?22 + 2^ 23 ) ■ 
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FIG. 1: The self-energy and vertex diagrams for calculating the renormalization constants. 



gb -> tW qq' -> tb qb -> q't 




FIG. 2: Tree level Feynman diagrams for associated production. 



gb -> tW 




FIG. 3: The Feynman diagrams for SUSY QCD corrections of associate production. 
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FIG. 4: The counter term diagrams for associated production. 




FIG. 5: The Feynman diagrams for s-channel and t-channel loop corrections, respectively. 
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FIG. 6: The K factors as functions of M g or M\n f° r PP ~^ at the LHC, where different curves 
correspond different tan/3, assuming: Mq = 150GeV, Aq = — 200GeV and fi > 0. 
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FIG. 7: The K factors as functions of or Mo for pp — > tW at the LHC, where different curves 
correspond different tan/3, assuming: Mu 2 = 40GeV, Aq = — 200GeV and \i > 0. 
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FIG. 11: The scale dependence of LO and SUSY QCD corrected cross sections of associated 
production at the LHC, we set the factorization and renormalization scales as /x/ = \i r = and 
/zq = mw + m t , assuming: tan/3 = 5, Mq = 150GeV, M]y 2 = 70GeV, Aq = — 200GeV and fj, > 0. 





FIG. 12: The K factors for the pp — > tb at the LHC. The variables are (a): Mg{ M^iq) and 
Mj- (Mq), respectively, assuming: tan/3 = 5, Aq = — 200GeV and fi > 0. 
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FIG. 13: The K factors for the pp — > tb at the LHC. The variables are (a): M^(Af 1 / 2 ) and 
M^(Mo), respectively, assuming: Aq = — 200GeV and fi > 0. 
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FIG. 14: The K factors as functions of (Mo) for the pp — > tb at the Tevatron. The graph shows 
different M\/2, assuming: tan/3 = 5, Aq = — 200GeV and fi > 0. 




FIG. 15: The K factors as functions of Mj(M^ 2 ) for the pp — > qt at the LHC, the graph shows 
different Mq, assuming: tan/3 = 5, Aq = — 200GeV and fj, > 0. 
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FIG. 16: The scale dependence of LO and SUSY QCD corrected cross sections of s-channel at the 
LHC, we set the factorization and renormalization scales as fif = /x r = /x, and /io = Tnt, assuming: 
tan/3 = 5, M = 150GeV, M 1/2 = 70GeV, A = -200GeV and n > 0. 
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FIG. 17: The scale dependence of LO and SUSY QCD corrected cross sections of t-channel at the 
LHC, we set the factorization and renormalization scales as [i$ = fi r = fi, and /j,q = mt, assuming: 
tan/? = 5, M = 150GeV, M 1/2 = 70GeV, A = -200GeV and fj, > 0. 



